Anomalously small wave tails in higher dimensions 



:3 
< 



Piotr Bizoii,^ Tadeusz Chmaj,^'^ and Andrzej Rostworowski^ 

' M. Smoluchowski Institute of Physics, Jagiellonian University, Krakow, Poland 
^H. Niewodniczanski Institute of Nuclear Physics, Polish Academy of Sciences, Krakow, Poland 
'"^Cracow University of Technology, Krakow, Poland 
(Dated: February 5, 2008) 

We consider the late-time tails of spherical waves propagating on even-dimensional Minkowski 
spacetime under the influence of a long range radial potential. We show that in six and higher even 
dimensions there exist exceptional potentials for which the tail has an anomalously small amplitude 
and fast decay. Along the way we clarify and amend some confounding arguments and statements 
in the literature of the subject. 

"iS": 

O . I. INTRODUCTION 

o ■ 

It is well known that sharp propagation of free waves along light cones in even-dimensional flat spacetimes, known as 
Huygens' property, is blurred by the presence of a potential. Physically, the spreading of waves inside the light cone is 
caused by the backscattering off the potential. If the potential falls off exponentially or faster at spatial infinity, then 
the backscattered waves decay exponentially in time, while the long range potentials with an algebraic fall-off give 
fT^ ■ rise to tails which decay polynomially in 1/t. The precise description of these tails is an important issue in scattering 
T-H ' theory. There are two main approaches to this problem in the literature. On the one hand, there are mathematical 
. results in the form of various decay estimates. These results are rigorous, however they rarely give optimal decay 
^ ' rates inside the light cone and provide very poor information about the amplitudes of tails. The notable exception is 
Q^, the work of Strauss and Tsutaya [l| (recently strengthened by Szpak where the optimal pointwise decay estimate 
^ for the tail was proved in four dimensions. Unfortunately, to the best of our knowledge, there is no analogous result 
in higher dimensions. 

On the other hand, there are non-rigorous results in the physics literature based on perturbation theory. The most 
^ ' complete work in this category was done by Ching et al. ^] who derived first-order approximations of the tails for 
^ I radial potentials. Although these results were originally formulated for partial waves in four dimensions, they can be 
' easily translated to spherical waves in higher dimensions. Ching et al. noticed that there are exceptional potentials for 
. which the first-order tail vanishes, however they did not pursue their analysis to the second order, apart from giving 
C ■ some dimensional arguments. The main purpose of this paper is to analyze the tails for such exceptional potentials 
'""^ [ in more detail. 

00 ■ One of the physical motivations behind our work stems from the fact that this kind of potentials arise in the 
[ study of linearized perturbations of higher even-dimensional Schwarzschild black holes. The behavior of tails on the 
C : Schwarzschild background in well known in four dimensions (see i, ^, %, i, 0, i), but not in higher even 
' dimensions (despite statements to the contrary in the literature [§]). Although our analysis is restricted to the flat 
J> . background, it sheds some light on the problem of tails on the black hole background because the properties of tails 
are to some extent independent of what happens in the central region. 
. , The rest of the paper is organized as follows. In section 2 we construct the iterative scheme for the perturbation 
' expansion of a spherically symmetric solution of the linear wave equation with a potential. This scheme is applied 
in section 3 to derive the first and second-order approximations of the tails for radial potentials which fall off as 
pure inverse-power at infinity. In section 4 we discuss the modifications caused by subleading terms in the potential. 
Section 5 contains numerical evidence confirming the analytic formulae from sections 3 and 4. Finally, in section 6 we 
give a heuristic argument to predict the behavior of tails outside Schwarzschild black holes in higher even dimensions. 
Technical details of most calculations are given in the appendix. 

Thr oug hout the paper we use the succinct notation and summation technics from the excellent book by Graham 
et al. [14| . In particular, we shall frequently use the following abbreviations 

x^:=x-{x-l) (x-(fc-l)), fc>0, (1) 

x°:=l, x'^ := X ■ (x + 1) (a; + (fc - 1)), fc>0. (2) 



X 



II. ITERATIVE SCHEME 

We consider the wave equation with a potential in even-dimensional Minkowski spacetime R'^^^ 

df(j) - A0 -f XV^ . (3) 
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The prefactor A is introduced for convenience - throughout the paper we assume that A is smah which allows us to 
use it as the perturbation parameter. The precise assumptions about the fall-off of the potential will be formulated 
below. We restrict attention to spherical symmetry, i.e., we assume that (f> — (j){t^r) and V — V{r). Then, equation 
([3]) becomes 

L(j) + \V{r)(j) = , C-.^d'i^dl- '^—^dr . (4) 

r 

We are interested in the late-time behavior of 4>{t, r) for smooth compactly supported (or exponentially localized) 
initial data. 

0(O,r) = /(r), dtc^{0,r) ^ g(r) . (5) 
To determine the asymptotic behavior of solutions we define the perturbative expansion (Born series) 

= ^ A"</)„ , (6) 

where (fiQ satisfies initial data ^ and all 0„ with n > have zero initial data. Substituting this expansion into 
equation ^ we get the iterative scheme 

£0^ = -F(/.„_i , 0_i=O, (7) 

which can be solved recursively. The zeroth-order solution is given by the general regular solution of the free radial 
wave equation which is a superposition of outgoing and ingoing waves [To| 

0o(t,r) = 0-*(t,r)+0g'''^(t,r), (8) 

where 

and u — t — r, v — t + r are the retarded and advanced times, respectively. Here and in the following, instead of d, 
we use the index I defined by d = 2Z -f 3 (remember that we consider only odd space dimensions d). Note that for 
compactly supported initial data the generating function a{x) can be chosen to have compact support as well (this 
condition determines a{x) uniquely). 

To solve equation ([7]) for the higher-order perturbations we use the Duhamel representation for the solution of the 
inhomogencous equation Ccf) = N{t,r) with zero initial data 

t t+r-T 

= / p'+'Pi{f^)N{T,p)dp, (10) 

\t~r-T\ 

where P;(^) are Legendre polynomials of degree I and p — {r^ + — {t — T)'^)/2rp (note that — 1 < /i < 1 within 
the integration range). This formula can be readily obtained by integrating out the angular variables in the standard 
formula (j) = C'^* * N where G^'^*{t,x) = (2^T'^~^^)~^&(t)S'^^^t'^ — |a;p) is the retarded Green's function of the wave 
operator in d 4- 1 dimensions (see, for example, [ll|). 

It is convenient to express (|10l) in terms of null coordinates rj = t — p and = t + p 



^(t^^) ^ TU-Tl-T / dU iC-vy^'m)N{v,Odv, (11) 



t-{-r t — r 
1 [ rlf [ 

2/+3j,/+l 

\t-r\ ~i 

where now ^ = (r2 + {£.-t){t- ri))/r{S, - 77). Using this representation we can rewrite the iterative scheme ([7]) in the 
integral form 

t+r t~r 

= / J {^^vy+'Piit^)V{piv,mn-iiv,Odri. (12) 

\t-r\ 

This "master" equation will be applied below to evaluate the first two iterates for a special class of potentials. It is 
natural to expect that for sufficiently small A these iterates provide good approximations of the true solution. 
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III. PURE INVERSE-POWER POTENTIALS AT INFINITY 

In this section we consider the simple case (below referred to as type I) when the potential is exactly V{r) = 
for r greater than some vq > 0. We assume that a > 2. The modifications caused by subleading corrections to the 
pure inverse-power decay of the potential will be discussed in section 4. 

A. Generic case 

We wish to evaluate the first iterate r) near timelike infinity, i.e, for r — const and t — > oo. Thanks to the fact 
that 4>Q{r/, ^) has compact support we may interchange the order of integration in (fT2|) and drop the advanced part of 
(/>o(?7, S,) to obtain 

oo t+r 

Mt,r) = -^^^^ J dr^ J {^-vy^'~"Pi{t^W{r,,OdC, (13) 

— OO t—r 

where we have substituted V = 2^{^ — r])~'^ . Plugging @ into (fT5|) . after a long calculation (see appendix A for the 
technical details), we get 



2 / V2 



[{t - 77)2 - r2 

n\{2l + 2n + iy. ^ 



0<n< L(q-2)/2J 



X 



2m 

r 



^ ' \m — nl (ay^ w _ jT 

n<m<n+l ^ ^ V 2 / 

Asymptotic expansion of (jl4p near timelike infinity yields the following first-order approximation of the tail 



0(t,r)«A0i(i,r)=A^^^ 



where 



A+ia + 2l)j +o(^ 



^(^'")^-(27TT)!!(^j (2) ' 

and 

+00 +00 

J a{7])dr], B= j a{iii)r]di^. (17) 

— oo — oo 

In general A^Q and the tail decays as t^"^'^', however there are nongeneric initial data for which A ^ and then the 
tail decays as ^-"-2/-!. particular this happens for time symmetric initial data for which a{x) is an odd fmiction. 

Remark 1. It is easy to check that if the function 4'{t,r) satisfies equation fl]), then the function — r^^^4) satisfies 
the radial wave equation for the lih multipole 

(5^ - dl + 1(1 + l)/r2)V + \Vir)iP = . (18) 

The late-time tails for this equation were studied by Ching et al. ^ who derived the formula equivalent to (llSp via 
the Fourier transform methods. 



(15) 
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B. Exceptional case 



It follows from (fT5|) that if a is an odd integer satisfying 3 < a < 2Z + 1, then r) vanishes identically due to 



factor in (|16p and there is no (polynomial) tail whatsoever in the first order. Thus, in order to compute the 

tail in this exceptional case we need to go the second order of the perturbation expansion. 
Using and proceeding as above we get the second iterate 



2" 



t+r 



(19) 



where (pY'* is the outgoing solution of the inhomogeneous equation 

£01 - -1/(^0 . 



(20) 



In general (j)i is a sum of the solution of the homogeneous equation and the particular solution of the inhomogeneous 
equation. The homogeneous part has the form ^ (with a different generating function than a, but still compactly 
supported), thus for the same reason as above it gives no contribution to the tail. The particular solution of the 
inhomogeneous equation ([20|) reads 



1 



2{a~ ly 



a+l 



E 

g=0 



(Z-a/2 + l/2)i 



29 (a/2)« 4>ti- 



(21) 



where (f>jLi^q denotes the solution of the homogeneous equation with d = 2(Z — 1 — g) + 3 and the same generating 
function a as in (pQ (see The formula (PT|) can be easily derived by the method of undetermined coefficients (we 
emphasize that this formula is valid only for odd a satisfying 3 < a < 2/ + 1). Substituting ([21]) into (I19p . after a 
long calculation (see appendix A for the technical details), we obtain the following asymptotic behavior near timelike 
infinity 



A + 2{a + l-l)j + o(^ 



(22) 



where the coefficients A and B are defined in (|T7|) and 

22{^+i-2) (2a -3)/ 5\ — , , 

^(''")=(2m)!!-kriy("-2j 

Here F stands for the generalized hypergeometric function 



f[ + "/2 - 1/2, a/2, 2a - 2, 1 
' ' a, a, a — Z — 1/2 



(23) 



F 



ai, 



(Xra 
bn 



k>0 "l' 



fofc fc! 



(24) 



TABLE I: The first few coeflicients D{1, a) 





3 


5 


7 


9 


1 


4 








2 


-8/5 


2240/3 






3 


96/35 


1792 


2523136/5 




4 


-64/7 


-17920/9 


16580608/5 


4638965760/7 



We remark that the behavior 0{\'^) t ^C+a i) gf ^^ji was conjectured before by Ching aL [3| on the basis 
of dimensional analysis. 
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IV. GENERAL POLYNOMIALLY DECAYING POTENTIALS 

In this section we analyze how the presence of subleading corrections to the pure inverse-power asymptotic behavior 
of the potential affects the results obtained in section 3. We restrict ourselves to the most interesting and common 
case (below referred to as type II) when near infinity 

y(-)-^{^ + ^)+o{^), 7>0. (25) 
If C(a, I) ^ 0, then the dominant behavior of the tail is of course the same as in ([T5|) : 

(t)it,r) ^ \ACil,a)t-^°'+^^K (26) 

However, in the exceptional case, when C(a, I) — 0, the situation is more delicate. As we showed above, in this case 
there is the second-order contribution to the tail given by (|22ll 

02(t, r) - AD{1, a) t-2(a+i-i) _ ^27) 

In contrast to the type I case where the first-order tail vanishes identically, in the type II case the subleading term in 
the potential produces the first-order contribution which is given by (|15p with a replaced by a -I- 7: 

01 {t, r)^ l3AC{l,a + -f) t-("+^+20 , (28) 

assuming that a -t- 7 is not an odd integer < d — 2 (otherwise one has to repeat the analysis for the next subleading 
term in the potential). 

Now, comparing the decay rates in (P7|) and we conclude that the leading asymptotics of the tail is given by 
the first-order term X(l>i{t, r) if 7 < a — 2 (we call it subtype Ila), but otherwise, i.e. for 7 > a — 2 (subtype lib), the 
second-order term \'^4>2{i, f) is dominant for t ^ 00. 

Remark 2. In the context of equation pS]) a formula analogous to ((28t was obtained by Hod who studied tails in the 
presence of subleading terms in the potential (see subgroup Illb in 12]). However, Hod's analysis, restricted to the 
first-order approximation, was inconclusive because, as we just have shown, without the second-order formula (j27|) 
one is not in position to make assertions about the dominant behavior of the tail. 

V. NUMERICS 

In order to verify the above analytic predictions we solved numerically the initial value problem (4-5) for various 
potentials and initial data. Our numerical algorithm is based on the method of lines with finite differencing in 
space and explicit fourth-order accurate Runge-Kutta time stepping. As was pointed out in [3], a reliable numerical 
computation of tails requires high-order finite-difference schemes, since otherwise the ghost potentials generated by 
discretization errors produce artificial tails which might mask the genuine behavior. The minimal order of spatial 
finite-difference operators depends on the fall-off of the potential - for the cases presented below the fourth-order 
accuracy was sufficient, but for the faster decaying potentials a higher-order accuracy is needed. To eliminate high- 
frequency numerical instabilities we added a small amount of Kreiss-Oliger artificial dissipation All computations were 
performed using quadruple precision which was essential in suppressing round-off errors at late times. 

The numerical results presented here were produced for initial data of the form 



0(0, r) = exp(-r2), 9*0(0, r) = exp(-r2) . 



(29) 
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As follows from ^ the generating function for these data is 

a{x) = 2"('+2' (1 - 2x) exp(-a;2) , hence A 
We considered the following potentials 



+ CXD 



tanh 



a+2 , 



Vir) 



tanh 



a+2 , 



1 



tanh''' r 
ri 



a{x)dx = -s/tt/S 

(type I) 
(type II) , 



1+2 



(30) 

(31a) 
(31b) 



for various values of a and 7. The regularizing factor tanh(r) introduces exponentially decaying corrections to the pure 
inverse-power behavior at infinity but such corrections do not affect the polynomial tails. The numerical verification 
of the formulae (|28p . and (|22p is shown in tables II and III. The observed decay rates agree perfectly with 

analytic predictions, while small errors in the amplitudes are due to (neglected) higher-order terms in the perturbation 
expansion. 

TABLE II: The generic case: numerical verification of the analytic formula (I15|l for the potential (31a) (A = 0.1) and initial 
data \29\ . Comparing the second column of this table (corresponding to a = 3.01) with the last column of table 111 one can 
see the discontinuity of the decay rate at a = 3 (for d = 5 and 7) . 





a = 2.5 


a = 3.01 


a = 4 


Theory 


Numerics 


Theory 


Numerics 


Theory 


Numerics 


Exponent 

d = 3 

Amplitude 


2.5 
-0.1253 


2.499 
-0.0881 


3.01 
-0.1785 


3.009 
-0.1518 


4 

-0.3545 


4.00002 
-0.3320 


Exponent 

d= 5 

Amplitude 


4.5 
0.0261 


4.501 
0.0235 


5.01 
-0.00089 


5.0101 
-0.00085 


6 

-0.2363 


5.9999 
-0.2318 


Exponent 

d= 7 

Amplitude 


6.5 
-0.0294 


6.501 
-0.0276 


7.01 
0.00089 


7.01 
0.00087 


8 

0.1418 


7.9999 
0.1404 



TABLE III: The exceptional case: comparison of analytic and numerical parameters of the tails for the potential (31b) (the 
first two columns) and (31a) (the third column) with a = 3, A = 0.1, and initial data ([29]). The analytic results are given by the 
formula p8|l for the subtype 11a potential, and by the formula (|27p for the type I and lib potentials. Note that although the 
dominant tails for the type I and the subtype lib potentials are theoretically the same, in the case lib there is an additional first 
order error due to the subdominant term 0{X)t~'''^^^"'^^^ which accounts for a slight difference in numerical accuracy between 
these two cases. 





7 = 1/2 (subtype Ila) 


7 = 1.75 (subtype lib) 


(type I) 


Theory 


Numerics 


Theory 


Numerics 


Theory 


Numerics 


Exponent 

d = 5 

Amplitude 


5.5 
-0.0731 


5.4993 
-0.0696 


6 

0.00886 


6.002 
0.00862 


6 

0.00886 


6.0000 
0.00843 


Exponent 

d = 7 

Amplitude 


7.5 
0.0603 


7.4998 
0.0579 


8 

-0.00177 


8.0003 
-0.00175 


8 

-0.00177 


7.9999 
-0.00172 


Exponent 

d = 9 

Amplitude 


9.5 
-0.1131 


9.4999 
-0.1115 


10 
0.00152 


9.9957 
0.00145 


10 
0.00152 


9.9997 
0.00149 
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VI. SCHWARZSCHILD BACKGROUND 



Consider the evolution of the massless scalar field outside the d + 1 dimensional Schwarzschild black hole 

-1 

„2 , ^2jo2 



ds^ = -^l--^jdt^+^l--^j dr^+r^d^U, (32) 

where dVi\_^ is the round metric on the unit sphere S'^^^ and d > 5 is odd. Here we use units in which the horizon 
radius is at r = 1. Introducing the tortoise coordinate x, defined by dr/dx = 1 — I/r'*^^, and decomposing the scalar 
field into multipoles, one obtains the following reduced wave equation for the jth multipole 13 1 

^,-^,,UM*.0. ^.(l--l,)(<?l±l^lRl±^,<ti)!). ,33, 

Note that p3p is the 1 + 1 dimensional wave equation on the whole axis — oo < a; < oo. For large positive x we have 

O -^], (34) 



d-3x'i-3 {2d - 5){d-3) x^d~5 



which implies that 



, , (2j + d - 3)(2j + d - 1) , , , , a b , ^ , 



with 



{d-l)j{j + d-2) {2d - 3)((d - 3)(d - 2nd - 1) ~ 4j(j + d - 2)(1 + d{d - 3))) 
" = d—3 ^= 4(2d-5)(d-3)^ • 

For large negative x (near the horizon) the potential is exponentially small, so one expects that the backscattering 
off the left edge of the potential can be neglected. If so, the decay rate (but not the amplitude!) should follow from 
the analysis of section 4. Comparing equation ((33|) for large positive x to equation (fT8|) with the potential ((25|) and 
using ([55)1 we find that I = j + {d — 3)/2 and the potential V is of the subtype Ila with a = d and 7 = d — 2. Thus, 
applying ([28|l we get the first-order tail 

^{t, x) - t-(2j+3d-5) _ (37) 

Remark 3. Late-time tails outside higher dimensional Schwarzschild black holes were studied in f^, however in the 
even-dimensional case the reasoning presented there is not correct, even though the result agrees with p7p . The reason 
is that the analysis of [9'] is based on the application of Ching et al. conjecture about the decay of the second-order 
tail i-(2'+2a-2) -^^hich for I = j + {d—3)/2 and a — d gives t^(2i+3'i-5). Unfortunately, this conjecture does not apply 
to the problem at hand. For j = this is evident because the leading term in V (proportional to a;"'') vanishes (since 
by ([36|) a = 0), while the subleading term (proportional to a;~^2d-2)) jg Qf generic type. For j > this follows from 
the fact that the potential is of the subtype Ila. Thus, for all j > the dominant (first-order) contribution to the 
tail comes from the subleading term in the potential. The agreement of the decay rate obtained in Q with ([37|) is 
accidental and due to the fact that the subdominant term in pS]) (not considered in [^) is on a borderline between 
subtypes Ila and lib. 

Admittedly, the handwaving argument leading to (|37p is far from satisfactory. Unfortunately, we have not been 
able to carry over the analysis from sections 2-4 in the case of equation ([55]) . There are two difficulties in this respect. 
First, in contrast to the spherical case, Huygens' principle is not valid for the free wave equation in 1 -f- 1 dimensions. 
Second, there is no natural small parameter in the problem. In the impressive tour de force work [7| Barack showed 
how to overcome these difficulties for a restricted class of initial data in four dimensions. It would be interesting to 
generalize Barack's approach to higher even-dimensional Schwarzschild spacetimes. 
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APPENDIX A 



Throughout the appendix we use the notation of 14] in which the square bracket around a logical expression returns 
a value 1 if the expression is true and a value if the expression is false: 



[condition] — 



1 if condition = true 
if condition ~ false 



In order to derive the asymptotic behavior of the iterates (jl3p and (|19p near timclikc infinity (fixed r and t — > co) we 
need to evaluate the following expression 



^ +00 i+r 



Tit.r- /3, L) 



k=o •' ■ /; 



-00 t~r 



where 



and 



- t){t - ij) + 



- 77) 



From ^ and (11311 we have 



and from and (PT|) we have 



Since 0(77) has compact support, it is advantageous to begin with integration by parts 

+00 +00 



a(?7) . 



— 00 —00 

For fi as defined in (jA3p and for any function g(fi) the following identity holds 

1=0 \ j 



hence 

+00 / \ t-\-r 



(A3) 



i{t,r)^T{t,r;a,l), (A4) 



^2 = 2(^_\w E - "/2 + 1/2)^ • ^(i, r; 2a - 1 + <z, / - 1 - g). (A5) 



nt^r;P,L) = -^ I dr^air,) ^ (-1)' T c.,.(/3 + L - 1)^1^ / ^^Ij^^p/^Ha^) • (A7) 



-00 u^j^ft,^j. \ ^ / t-r 
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1. First-order approximation 

To evaluate the first iterate (j>i{t,r) we apply the formula (jA13|l with (3 — a and L = I. Then 

where we shifted the summation index m j + n — m. Next, we interchange the order of summation according to 

[0<j][j<l][rn-n<j][j<2m + l] 
^ [-^ < "I < 0][0 < j < ; + 2m] + [0 < m < 7i][0 < j < + [n < m < I + n][m - n < j < I] , 

and convert the sum over j into a generahzed hypergeometric function [14]. Defining 



a-2\ll + 2n + j\ l j \ I a~2 + l~j 
\ j \ j + n-m I \ l-j + 2m 



we see that io 7^ iff n = to, thus the sums for [— ^ < m < 0] and [0 < m < n] do not contribute to (|A15P and we are 
left with 



n+l / \ 2m l+n—m 

m—n ^ ^ j—0 



a — 2 \ll + n + m + j 
I + n — m — j I \ j + m — n 



j + m — n \ I a — 2 + l + n — m — j 
j \ l + n + m- j 



(A16) 



where we shifted the summation index j ^ j + m — n. Defining 

i, = (-1)^ 



a — 2 \ I I + n + m + j \ I j + m — n \ I a — 2 + l + n~m~j 
I + n — m — j I \ j + m — n j \ j j \ I + n + m — j 



we see that 



and 



hence 



~ _ (g- 2)'+"-" {l + n + my. {a ~ 2 + I + n - to)'+"+'" 
(Z + n — to)! (to — n)!(/ + 2n)! {l + n + my. 

Ij+i _ {j - {l + n- m)){j ~{l + n + m)){j + {l + n + m+l)) 

1~ ^ (j + i{a-l)-{l + n- m))){j + {-{a - 2) - [l + n - to)))(j + 1) ' 

(q - 2) '+"-"' {a- 2 + l + n- to) '+»+"' 

—n ^ ' ^ 



r\ ) (Z + n — to)! (to — n)!(Z + 2n)! 

— (l + n — to), —(l + n + to), (I + n + m+X) 
(a — X) — (I + n — m), —{a — 2) — (l + n — to) 

/ X2m ^M+n-m (q, _ 2 + / + n - to)'+"+"' 



E[ ^ ) gl+2(i+n-m)^ jO^- ^J- 
\t-nj il + n- 



T] J {I + n — my. {m — ny.(l + 2ny. 

T{-{a -2)- {l + n- TO))r((a - 1) - (Z + n - to)) 

r (-^ + ™) r (-^ -{1 + n)) r (f + to,) r -{i + n)) ' 



(A17) 
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where in the last equation we used the identity 



F 



a+l, -a, {b + c- l)/2 
b, c 



1 = 22-(''+-)7r 



mm 



r(^)r(^)r(i±|±^)ri 



Substituting 



and 



(a _ 2)^±I^r((a - 1) - (/ + n - m)) = r(a - 1), 



{a~2 + l + n- my-±^^r{-{a -~ 2) - {I + n - m)) = (-l)'+"+™r(-a + 2 + 2m) 



into (|A17[) we get 



n+l 



^ \2ni ^_-[^y+n+m2l+2(i+n-m)^ 



r(a- l)r(-a + 2 + 2m) 



r^y (/ + n - m)!(m - n)!(Z + 2n)! r (f + to) T (-^ + m) T (-^^3^ - + n)) T - {I + n)) 
The last equation can be still simplified due to the identity 



r(a- l)r(-a + 2) 



(-1)' /a-3\-/a 



We have 



r(f)r(-^)r(-^-/)r(^ifi-0 2n \ 2 



T{-a + 2 + 2m) = (-a + 2)2"T(-a + 2) 



(f) 



a — 3 



a — 3\ ^/ a — 3 



a-2 
2~ 

a- 1 



Z — n 



I - n 



2 



a-2 



-l-l)- 



r ('I + m 



(^_-;-if 

(2) ^(2)' 



and 



/ a^V" V 2 / 



so finally 



n+l 



2m 



^ -j^ ^n+m 22 (/+n) 



77/ (/ + n - to)!(to - + 2n)! V 2 
irf-^-/-n-f^i^-?-l)- 



(f) 



(f)'" 

Plugging (|A19p into (|A13p with (3 = a and L = I we get the expression ([14]). 



(A18) 



(A19) 
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2. Second-order approximation 



The calculation in the second order {f3 — 2a ~ 1 + q and L — I — 1 — q) is only a slight modification of what we 
have already done in the first order. Following the same steps which led us from (jAlSp to (jAlSp we get 

"+-^ / ^ \2ni ^_^y+n+rn2l+2{L+n-m,)^ 



-^^_^\t — rj J (L + n — m)!(m — n)!(Z + 2n)! 

r(/3 - l)r(-/3 + 2 + / - L + 2to) 

The last equation can be simplified due to the identity 

r(/3 - i)r(-/3 + 2 + /-L) 



L / a , 1 T \L 



{-ly fp-3-{l-L)\- fl3 + l-L 



(/3-2)^, (A21) 



2tt 

which for L = I reduces to (jAlSp . We have 

T{-(3 + 2 + l-L + 2m) = {-/S + 2 + I - Lf^T{-l3 + 2 + I - L) 
P( ^~^ + ''^^+rr?j-(^ /3-3, l-L\^^f /3 - 3 , l-L 

r 



(3-2 l + L 
1 ^ :^ n 



( 0-2 _ 

;3-2 _ J+L ^ \31 ' 
2 2 



/3~ 1 ; + L \ ^ 2- 



r ; ^ n 



X n ■) 

/3-1 _ /+L _ 1 
2 2 



and 



hence 



2 



2 2 



- ^ {L + n - my.{m, - ny.{l + 2ny. 



(A20) 



(A22) 



2 

Plugging (|A22p into (|A13|) we get 



( p .l-L 
\^2 2 



T{t,r- 2a-l + qJ-l-q) = (-1)' « - (a - 2 + (2a - 3)1+« 



2^ 

— oo 

^^+^-1-'? / 7 1 ^ \ / 1 / \ 2m 



n=0 



' ■ ^' n!(2/ + 2n + l)! ^ , ^ ~, ^ 2 ^ ' V ^„ ' " I m - n / + V^-^. 



r 



13 



Substituting this into (|A5|) and expanding in 1/t we have 

2 22"+2'-2 2 



2(a-l) 4(2/ + 1)!! t2a+2/-2 



A + 2(a + /-l)| 



9=0 



l-l-q 



"(q- 2 + 0^^^ (2a- 3)1+* 



with and B defined in p7p . Converting the sum over g into the generahzed hypergeometric function we get 



(A24) 
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